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We consider the relation between relaxation time and the largest Lyapunov exponent in a system 
of two coupled oscillators, one of them being harmonic. It has been found that in a rather broad 
, region of parameter space, contrary to the common expectation, both Lyapunov exponent and 

f~>) ■ relaxation time increase as a function of the total energy. This effect is attributed to the fact that 

' above a critical value of the total energy, although the Lyapunov exponent increases, KAM tori 

appear and the chaotic fraction of phase space decreases. We examine the required conditions and 
demonstrate the key role of the dispersion relation for this behaviour to occur. This study is useful, 
among others, in the understanding of the damping of nuclear giant resonances. 

PACS numbers: 05.45.Pq, 24.60.Lz, 05.45.Xt, 24.30.Cz 
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I. INTRODUCTION 
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. It is well known that time-independent Hamiltonian systems with at least two degrees of freedom can exhibit chaotic 
behaviour 0,0- In most cases the dynamics is mixed, i.e. regular and chaotic regions coexist in phase space. In 
general, the largest Lyapunov exponent (Ai) is the most common measure of the chaoticity of the system. It appears 
that, at least in completely chaotic systems, increasing Ai, the relaxation time, defined as the time required for an 
observable to reach its equilibrium value, decreases 0,G1- However, in mixed systems the relationship between Ai and 
relaxation time is not well established yet. 

Relaxation time is an important tool in many body problems [E IS 0> HI 0> Ell EH EH El El E3 ■ Let us suppose 
that a many body system is collectively excited by an external perturbation which for example may be a collision 
with a particle-projectile or a laser pulse. The energy of the collective excitation is gradually absorbed by the internal 
dynamics of the system and the time needed for this absorption to occur is the relaxation time of the excitation. 
' Such relaxation processes can be encountered for example in the field of nuclear physics and in particular in the 
damping of nuclear giant resonances 0, 0, 0|. In the latter the system has mixed dynamics. The mechanism 
responsible for the dissipation of nuclear collective energy is still an open problem. Especially, the effects of chaotic 
dynamics on the damping of nuclear excitations is still not completely understood. One body dissipation |19| | refers 
to processes involving a single nucleon (with mean free path of the order of the size of the nucleus) interacting 
with the collective nuclear potential. We have employed a simple classical model based on this independent particle 
approach for the dynamics of the nucleons. This model can be used for the interpretation of the isoscalar monopole 
resonance |2(| (the breathing mode) decay 0, ITU l2ll |2^| . It consists of a harmonic oscillator describing the collective 
excitation (to a first approximation) coupled with a nonlinear (Woods-Saxon) oscillator, which, in the independent 
particle approach, represents the motion of each nucleon in the average potential created by all nucleons. The model 



corresponds to a time-independent Hamiltonian system with 2 degrees of freedom. The detailed investigation of the 
specific giant resonances through this simple classical model will be given elsewhere [23|. Here, we will use this model 
in order to elucidate the relationship of relaxation dynamics with the structure of the chaotic component of the phase 
space of the system. We should add that an appropriate classical model has also been used to study the damping of 
giant dipole resonance in Ref . [24J . While the intensity of that chaotic component is usually measured by the largest 
Lyapunov exponent as mentioned above, its extent is quantified by the fraction of phase-space occupied by chaotic 
orbits. The aim of this paper is to explore the interrelation between the relaxation time of a harmonic oscillator and 
Ai, as well as how this is affected by changes in the chaotic fraction of phase space. 

More specifically, as the considered model posseses mixed phase space, we have found that, contrary to a completely 
chaotic system, a broad range of parameter space exists where Ai increases as the total energy increases and also 
both Ai and the relaxation time of an ensemble of chaotic orbits increase. We demostrate that this relationship is 
affected by the form of the dispersion relation lu(E) of the nonlinear oscillator. In particular, the frequency ui(E) 
must either be a strictly increasing function of the energy or have a maximum (zero dispersion oscillator |25j) for 
the above mentioned behaviour to occur. In these cases, the chaotic fraction of phase space is shrinking as the total 
energy increases, inducing a delay in the relaxation of the dynamics. Finally, the dependence of the results on the 
form of the coupling as well as on the particular choice of the potential of the non-harmonic oscillator is explored. 
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The paper consists of four further sections. In Sec. II we describe our model system. In Sec. Ill we present our 
findings on the relation between relaxation time and Lyapunov exponent for a particular case. In Sec. IV we discuss 
the dependence on the system parameters of the effects presented in Sec. III. In Sec. V we summarize our main 
results and give future prospects. 

II. DESCRIPTION OF THE MODEL SYSTEM 

Our model system consists of a particle moving in a Woods-Saxon well of finite depth V coupled with a harmonic 
oscillator moving with frequency loho ■ The Hamiltonian of the system is 



where the potential V(r, R) is given by 



V ^ ^ = (r V ° R +U R KT\ + \ M ^ {R R0)2 - (2) 

i + cx P r~ K ° +fc a (Ko '- R) j 2 

A contour plot of the potential for a specific choice of the parameters, namely a = 0.5/m, b = 0.2, i? = 7.69/to, 
V a = 48MeV, M/m = 1 and uj = 13.73MeV/h, is shown in Fig.l. 

The choice of this Hamiltonian is motivated by the damping of a density collective mode in nuclei. Specifically, this 
mode is represented by a harmonic oscillator of mass M which interacts with the nucleons, represented as independent 
particles of mass m moving in an average potential which has the form of a Woods-Saxon well. The strength of the 
coupling between the particle and the harmonic oscillator is controlled by the parameter b. In the limit 6^0 the 
oscillators become uncoupled and thus the system becomes integrable. As b increases the system becomes increasingly 
chaotic. Hamilton's equations of motion for the system are 

^(z=gzg) „ 



/ r-R o +b(R -R) \ 

PR = \ " ^72 - (R Ro) . (6) 



a (l + cw (r-R o+b (R -R) ^y 



These equations are solved with the use of an adaptive fourth order Runge-Kutta algorithm. The total energy E was 
conserved with relative error AE/E < 10~ 5 . Rescaling the differential equations of motion we find that the relevant 
parameters of the system are M/m, Ro/a, Vb/(ma 2 Wp) and b. 

III. RELAXATION TIME, LYAPUNOV EXPONENT AND CHAOTIC FRACTION OF PHASE SPACE 

For the analysis of the system dynamics, we start with a microcanonical ensemble of initial conditions (5000 in 
most of the cases) in a rectangular region of (r,E p ) plane, where E p is the particle energy defined as 

vl v 
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FIG. 1: Contour plot of the potential V(r,R) of the system for a = 0.5 fm, b = 0.2, Ro = 7.6 fm, Vo = 48MeV and 
ujo = 13.73MeV/h (see Eq. 2). 



For each such ensemble, the total energy E = E p + Eho is constant, where Eho is the energy of the harmonic 
oscillator given by 

E HO = l -Mu J l{R~R f + ^ } . (8) 

All initial conditions we use belong to the chaotic region of phase space. We evolve them forward in time and at each 
time instant t we calculate the mean value < R(t) > of R. It is expected that < R(t) > will approach the equilibrium 
value R eq (which is very close to Rq) almost exponentially since all orbits in the ensemble are chaotic. The 
corresponding time scale defines the relaxation time. For the values of the parameters chosen, the time dependence 
of < R(t) > —R eq can be fitted very well by a function of the form R — R eq — Ae~ 7t cos (cut). We calculate the 
half-life Ti/ 2 of the initial amplitude A as a function of the total energy of the system. T^/ 2 and 7 are related by 
T1/2 — In 2/7. The structure of the phase space allows us, as it will become clear later, to use the same rectangular 
region of initial conditions for all the considered values of the total energy. However, we have found that there is no 
dependence of the relaxation time on the location of the distribution of the initial conditions in phase space, provided 
that this distribution lies entirely in the chaotic regime. 

To quantify the degree of chaoticity of these orbits, we calculate for each one the Lyapunov exponent Ai by averaging 
the growth rates of small perturbations along them |]J,|2|. It is found numerically that, for sufficiently large integration 
times (of the order of the relaxation time), the average value < Ai > calculated in this way is independent of the 
location of the rectangle of the initial conditions in the chaotic region. We determine the average value < Ai > over 
the ensemble of initial conditions as a function of the total energy of the system. 

In addition to Ai we calculate the percentage q of chaotic phase space in the (R,pr) plane as a function of the 
total energy of the system [2(||23]. In order to do this, we start with an ensemble of a large number of chaotic initial 
conditions on the (R,pa) section and we integrate them forward in time until they cover the entire accessible chaotic 
region. We then divide the (R,pu) section in N — 10 4 cells and count the number of occupied cells N occ , i.e. the 
number of cells in which there is at least one orbit point. We estimate the percentage of chaotic phase space as 

<z = ^pioo%. 

In the study presented in this section we choose the following values for the parameters of the system: Rq = 7. 69 fm, 
Vq = 48MeV and uiq — 13.73MeV/h. These parameters are relevant for the decay of the isoscalar giant monopole 
resonance of the nucleus 208 Pb H, |2(], HU . The strength of the coupling has been chosen equal to b = 0.08. We 
also choose M/m = 1 and a = 0.05/m. The results for < Ai >, T1/2, and q as a function of the total energy of 
the system are shown in Fig. 2. It is clearly seen that although the value of < Ai > is a monotonically increasing 
function of the total energy E (see for instance Q), the relaxation time, above a critical value of the total energy, stops 
to decrease and starts increasing contrary to the common expectation from fully chaotic systems. We also observe 
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FIG. 2: (a) < Ai >, (b) T1/2 and (c) q as a function of the total energy E for M/m — 1, Ro — 7.69/m, a — 0.05/m, 
V = 48MeV, lo = 13.73MeV/h and b = 0.08. 



FIG. 3: (R,pb) Poincare sections for (a) E = -35MeV, (b) E = -25MeV, (c) E = -15MeV and (d) E = —5MeV. The 
other parameters are the same as those of Fig. 1. 



that above this critical energy, the percentage q of chaoticity in the (R,pr) section is decreasing. Motivated by this 
observation we consider in more detail the Poincare surfaces of section for different values of the total energy of the 
system. In Fig. 3, (R,pr) sections (for r = Rq/2 and p r > 0) for four values of the total energy (E = — 35MeV, 
—25MeV, —15MeV and — 5MeV) are shown. In Fig. 4 we plot (r,E p ) sections (for R = R , pr > 0) for the same 
values of the total energy. 

The regular region of Figs. 3(c) and (d) contains invariant tori and corresponds to the upper invariant curves of 
Figs. 4(c) and (d) respectively. These curves are close to straight lines and correspond to almost constant particle 
energy and hence almost constant total energy. They are therefore Kolmogorov- Arnold- Moser (KAM) tori pj, 
which are remnants of the uncoupled integrable system. By inspecting these Poincare sections, we observe that the 
appearance of KAM tori leads to the observed relative decrease of chaos in phase space. The onset of this decrease 
coincides with a change in the monotonicity of Ti/ 2 as a function of the energy. Once the KAM tori appear, they are 
not destroyed as the energy increases. We conjecture that such an appearance of KAM tori counteracts the increase 
of the Lyapunov exponent and leads to a change in the monotonicity of the relaxation time. But what causes this 



FIG. 4: (r, E p ) Poincare sections for (a) E = —35MeV, (b) E = —25MeV, (c) E = —15MeV and (d) E = —5MeV. The other 
parameters are the same as those of Fig. 1. 
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FIG. 5: (a) The Woods-Saxon potential for Vb = 4SMeV and four values of a, namely 0.05/m, 0.5/m, 1.5/m, and 3.1/m, 
(b) The ujp(Ep) curve of the Woods-Saxon potential for the same values of the a. The straight line shows the frequency of the 
harmonic oscillator. 



reappearance of KAM tori in the system phase space for such large values of the energy? How is it connected to the 
basic properties of the non-harmonic oscillator ? To answer this question, we investigate the form of the dispersion 
relation of the non-harmonic oscillator. This relation, which is sensitive to the parameters of the system and controls 
the appearance of KAM tori, will be studied in the following section. 



IV. THE KEY ROLE OF DISPERSION RELATION AND THE INFLUENCE OF THE COUPLING 

STRENGTH 

The diffuseness parameter a is of crucial importance for the appearance of the effect described in the previous 
section. In the Woods-Saxon potential Vws( r ) = —Vo/ (l + cxp ( r ~ a R ° )) the frequency of the motion of the particle 
uj p depends on its energy E p , whereas the frequency of the harmonic oscillator u>ho is independent of the energy. The 
frequency of the motion of a particle in a potential V(r) is given by 

^'-'(/ v»(*%r)) ) ' <9) 

where ry and n are the right and left turning points at a given energy E p , i.e. the roots of the equation E p = V(r). 
The form of the u p (E p ) curve of the Woods-Saxon potential depends strongly on a as can be seen from Fig. 5, where 
V(r) and uj p (E p ) curves for four values of a are shown. At the limit a — > the potential tends to the square well and 
therefore the w p {E p ) curve has the form ui p cx ^f~E~p. For large values of a, uj p (E p ) is decreasing and for intermediate 
values it has a maximum. 

For a fixed value of a (i.e. a = 0.05/m) for which the observed behaviour of relaxation time occurs, the uu p {E p ) 
curve has two intersection points with the corresponding curve of the harmonic oscillator. The leftmost intersection 
point is a oj p : u>ho = 1:1 resonance. This resonance is apparent in Figs. 3 and 4. The rightmost intersection point 
is a much less apparent 1 : 1 resonance located close to the edge of the potential well. For small values of b, the 
phase-space is dominated by KAM tori. Resonances of small width also exist. As the coupling strength b increases, 
the width of these resonances also increases. They therefore overlap with higher-order resonances located close to 
them. This overlap leads to the appearance of a chaotic region. The area of this chaotic region increases with 6. 
Higher-order resonances (1 : 2, 1 : 3, ...) appear in the low-energy region and their density increases as the energy 
becomes smaller. Therefore, as b increases chaos emerges first in the low-energy region of the (r, E p ) section. This 
transition is illustrated in Fig. 6, where (r, E p ) sections for a = 0.05 fm and several values of b are shown. 

For fixed values of b and a, above some particle energy threshold only KAM tori exist. As the energy increases, the 
existing KAM tori are not destroyed and new KAM tori are added to the phase space in the region of large particle 
energies. This happens provided the coupling is not strongly chaotic. Therefore, as the energy increases, the relative 
area of the chaotic phase space decreases. This can be seen from Fig. 4. Increasing b makes the system more chaotic, 
i.e. both the Lyapunov exponent and the relative area of the chaotic phase space increase (see Fig. 6). The inverse 
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FIG. 6: (r,E p ) Poincare sections for E = -5MeV, a = 0.05/m and (a) b = 0.02, (b) b = 0.04, (c) b = 0.1, (d) b = 0.15, (e) 
b = 0.2 and (f) b — 0.3. The other parameters are the same as those of Fig. 1. 
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FIG. 7: (a) < Ai >, (b) Ti/ 2 and (c) q as a function of the total energy for a = 0.05/m and b = 0.08, b = 0.1 and b = 0.2. The 
other parameters are the same as those of Fig. 1. 



behaviour of relaxation time and Lyapunov exponent appears as long as a sufficiently wide layer of KAM tori exists 
at high energies. This can be seen from Fig. 7, where the Lyapunov exponent increases, the percentage of relative 
chaotic area decreases and the relaxation time, above a critical value of the total energy, either increases or remains 
almost constant. 

Increasing b further, will lead to a destruction of all KAM tori and thus to the usual behaviour of the relaxation 
time: < Ai > is an increasing function of the energy, q is also increasing due to the decrease of the area of the 
regular island around the 1 : 1 resonance and Ti/ 2 is decreasing. As b increases further (b = 3 and above), although 
the portion of the chaotic phase space grows, some stickiness effects may occur which also leads to an increase of 
relaxation time p7j . 

Increasing a leads to a decrease of the maximum value of the frequency of the Woods-Saxon potential and to a shift 
of its location towards lower energies. In Fig. 8, (r, E p ) surfaces of section for b = 0.2 and four values of a, namely 
0.05/m, 0.5/m, 1.5 fm and 3.1/m, are shown. These correspond to the three different qualitative forms of the uj p (E p ) 
curve (see Fig. 5). For a = 0.05/to and a = 0.5/m the curve has a maximum and two intersection points with the 
corresponding curve of the harmonic oscillator. The leftmost of these points corresponds to the 1 : 1 resonance which 
is located around E p = —ZQMeV and can be seen in Figs. 8 (a) and (b). For a — 1.5/m, uj p (E p ) has a maximum 
but no intersection points with the corresponding curve of the harmonic oscillator. In this case, as can be seen from 
Fig. 8 (b), chaos emerges first in the region of high particle energies. This is due to the fact that the u> p {E p ) curve 
is not symmetric with respect to its maximum: its slope is larger for energies above the maximum. Therefore, the 
density of the higher-order resonances (1 : 2, 1 : 3, . . .) is larger for these energies, and these resonances overlap, first 
creating a chaotic region at large energies. For a = 3.1, lo p (E p ) is decreasing and has one intersection point with the 
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FIG. 8: (r, E p ) Poincare sections for b = 0.2 and (a) a — 0.05/m, (b) a — 0.5 fm, (c) a — 1.5 fm and (d) a = 3.1 fm. The other 
parameters are the same as those of Fig. 1. 



corresponding curve of the harmonic oscillator. In this case, the 1:2,1:3,... resonances appear at higher energies 
and they become denser near the edge of the potential well, where the slope of the uj p (E p ) curve becomes large. 

The effect where Lyapunov exponent and relaxation time both increase can also be observed in simpler systems, 
provided that the dispersion relation lo(E) of the non- harmonic oscillator has an extended increasing part. Examples 
of such systems are oscillators with potentials ~ x 2n (n — 2, 3, . . .) coupled with harmonic oscillators. However, the 
generality of this effect is limited by the form of the coupling term. Actually it has been found that coupling terms 
exhibiting strongly chaotic behaviour (such as ~ x 2 y 2 ), destroy the above mentioned effect and the usual behaviour 
of relaxation times is recovered. 



V. CONCLUSIONS 



Motivated by questions of the damping mechanisms of collective motion in nuclei we have studied the behaviour 
of the relaxation time of an ensemble of chaotic orbits as a function of the total energy in a system consisting of 
two oscillators, one of them being harmonic. It has been found that, in a rather broad range of the parameters, the 
relaxation time can increase or remain constant although the Lyapunov exponent increases. We have attributed this 
behaviour to the appearance of KAM tori, which limit the phase space available to chaotic orbits. For this behaviour 
to occur, the dispersion relation of the non- harmonic oscillator must either be strictly increasing or have a maximum. 
The dependence of this behaviour on the parameters, as well as on the coupling has been investigated. In addition, 
it has been found that this behaviour can occur in very simple systems, provided that the coupling is not strongly 
chaotic. The model used in this work can be extended to study the direct decay of the isoscalar giant monopole 
resonances, i.e. the decay where nucleons can escape from the nucleus. A systematic study of this decay in several 
nuclei is presently under investigation p3| . 
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